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On the basis of a nine-parameter expanding source model that includes special 
C \ relativity, quantum statistics, resonance decays, and freeze-out on a realistic 

hypersurface in spacetime, we analyze in detail invariant vr^, vr^, , and 

■ K~ one-particle multiplicity distributions and vr"*" and two-particle cor- 
relations in nearly central collisions of Si + Au at p\a\^/A = 14.6 GeV/c. By 
considering separately the one-particle data and the correlation data, we find 

\ that the central baryon density, nuclear temperature, transverse collective 

^ \ velocity, longitudinal collective velocity, and source velocity are determined 

■ primarily by one-particle multiplicity distributions and that the transverse 

■ radius, longitudinal proper time, width in proper time, and pion incoherence 
00 ■ fraction are determined primarily by two-particle correlations. By consid- 



ering separately the pion data and the kaon data, we find that although 
the pion freeze-out occurs somewhat later than the kaon freeze-out, the 99% 
confidence-level error bars associated with the two freeze-outs overlap. By 
constraining the transverse freeze-out to the same source time for all points 



O 

P3 . with the same longitudinal position and by allowing a more flexible freeze-out 

in the longitudinal direction, we find that the precise shape of the freeze- 
^ \ out hypersurface is relatively unimportant. By regarding the pion and kaon 

' one-particle data to be unnormalized, we find that the nuclear temperature 

^ increases slightly, but that its uncertainty increases substantially. By includ- 

ing proton one-particle data (which are contaminated by spectator protons), 
we find that the nuclear temperature increases slightly. These detailed stud- 
ies confirm our earlier conclusion based on the simultaneous consideration 
of the pion and kaon one-particle and correlation data that the freeze-out 
temperature is less than 100 MeV and that both the longitudinal and trans- 
verse collective velocities — which are anti-correlated with the temperature — 
are substantial. We also discuss the flaws in several previous analyses that 
yielded a much higher freeze-out temperature of approximately 140 MeV for 
both this reaction and other reactions involving heavier projectiles and/or 
higher bombarding energies. 
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I. INTRODUCTION 



Experimentalists around the world are vigorously searching for the quark-gluon plasma — 
a predicted new phase of nuclear matter where quarks roam almost freely throughout the 
medium instead of being confined to individual nucleons . Such a plasma is believed to 
have existed in the first 10 /is of the universe during the big bang and could be produced in 
the laboratory during the little bang in a relativistic heavy-ion collision. 

When nuclei collide head-on at relativistic speeds, the nuclear matter is initially com- 
pressed and excited from normal nuclear density and zero temperature to some maximum 
values — during which pions, kaons, and other particles are produced — and then expands, 
with a decrease in density and temperature. The early stages of the process are often 
treated in terms of nuclear fluid dynamics, but at some late stage the expanding matter 
freezes out into a collection of noninteracting hadrons. 

Measurements of the invariant one-particle multiplicity distributions and two-particle 
correlations for the pions, kaons, and other particles that are produced sample the density, 
temperature, collective velocity, size, and other properties of the system during this freeze- 
out. The use of two-particle correlations to extract size information was pioneered in the 
1950s by BrownQ and Twiss 0, who used two-photon correlations to measure the size of 
stars, and by Goldhaber et al. IQ, who used two-pion correlations to measure the size of 
the interaction region in antiproton annihilation. The hope is that a sharp discontinuity in 
the value of one or more of the extracted freeze-out properties as a function of bombarding 
energy and/or size of the colliding nuclei could signal the formation of a quark-gluon plasma 
or other new physics. 

For the extraction of the freeze-out properties from experimental measurements of in- 
variant one-particle multiplicity distributions and two-particle correlations, a nine-parameter 
expanding source model that includes special relativity, quantum statistics, resonance de- 
cays, and freeze-out on a realistic hypersurface in spacetime was developed in Refs. P,^. 
The application of this model to central collisions of Si -t- Au at piabM = 14.6 GeV/c 
1^,11] led to the conclusion that the freeze-out temperature is less than 100 MeV and that 



both the longitudinal and transverse collective velocities — which are anti-correlated with the 
temperature — are substantial. Similar conclusions concerning a low freeze-out temperature 
have also been reached in Refs. |T^JT^. However, other analyses [^^-|TE[ have yielded a much 
higher freeze-out temperature of approximately 140 MeV for both this reaction and other 
reactions involving heavier projectiles and/or higher bombarding energies. 

Because of the importance of resolving this significant discrepancy, we perform here a 
more detailed analysis of exactly the same data |ll0|,p!T| that were considered in Refs. [BIH- 



^Like many individuals throughout history, R. H. Brown preferred his middle name to his first 
name. Therefore, instead of using Robert H. Brown as author on his many articles and books, he 
used R. Hanbury Brown instead. This practice inevitably led to the citation of Brown as Hanbury 
Brown and even Hanbury-Brown. Although this error is widespread throughout the relativistic 
heavy-ion community, he is correctly listed under the last name Brown in most astronomy texts, 
biographical reference books, and encyclopedias. 
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After briefly reviewing the source model in Sec. || and the results of a simultaneous con- 
sideration of the one-particle data and correlation data in Sec. Ill A , we examine various 



subsets of the data and determine the effect of alternative assumptions on our results. In par- 



ticular, we consider separately the one-particle data and the correlation data in Sec. [Ill B 
and consider separately the pion data and the kaon data in Sec. [Ill C 



In Sec. nm 



we 

constrain the transverse freeze-out to the same source time for all points with the same 
longitudinal position and allow a more generalized freeze-out in the longitudinal direction, 
and in Sec. [Ill E[ we regard the pion and kaon one-particle data to be unnormalized and also 
include proton one-particle data (which are contaminated by spectator protons). These de- 
tailed studies provide a useful background for our discussion in Sec. |rV| of the flaws in several 
previous analyses that led to an anomalously high freeze-out temperature. Our summary 
and conclusion are given in Sec. |V]. 



II. NINE-PARAMETER EXPANDING SOURCE MODEL 



The expanding source model introduced in Refs. describes invariant one-particle 
multiplicity distributions and two-particle correlations in nearly central relativistic heavy- 
ion collisions in terms of nine parameters, which are necessary and sufficient to characterize 
the gross properties of the source during its freeze-out from a nuclear fluid into a collection 
of noninteracting, free-streaming hadrons. The values of these nine parameters, along with 
their uncertainties at 99% confldence limits, are determined by minimizing for the types 
of data considered. Several additional physically relevant quantities, along with their un- 
certainties at 99% confidence limits, can then be directly calculated. The nine independent 
source freeze-out properties that we consider here are the central baryon density n, nuclear 
temperature T, transverse collective velocity ft, longitudinal collective velocity ve, source 
velocity Vs, transverse radius Rt, longitudinal proper time Tf, width in proper time At, and 
pion incoherence fraction A,r- 

For a particular type of particle, the invariant one-particle multiplicity distribution and 
two-particle correlation function are calculated in terms of a Wigner distribution function, 
which is the phase-space density on the freeze-out hypersurface, giving the probability of 
producing a particle at spacetime point x with four-momentum p. It includes both a direct 
term |]19[ and a term corresponding to 10 resonance decays |^, namely the decay of meson 



resonances with masses below 900 MeV and strongly decaying baryon resonances with masses 
below 1410 MeV. 

The direct part of the Wigner distribution function for a particular type of particle is 
given by g|TP| 



Q ( ) = '^£^ p-n{x) , . 

dir[x,p) (27r)3 exp{[j9-t;(x) -Ai(x)]/T(x)}Tl ' 

with the minus sign applying to bosons and the plus sign to fermions. The quantity J is 
the spin of the particle, v{x) is the collective four- velocity, T(x) is the nuclear temperature, 
and fi{x) is the chemical potential for this type of particle. The four- vector n{x), with 
components 
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n,{x)= / d\{x')5^^\x-x') , (2) 



E 



gives the normal-pointing freeze-out hypersurface elements. The subscript S on the integral 
denotes the limits to the hypersurface for a finite-sized system. Because we are considering 
nearly central collisions, we assume axial symmetry and work in cylindrical coordinates in 
the source frame, with longitudinal distance denoted by z, transverse distance denoted by 
p, and time denoted by t. Throughout the paper we use units in which h = c = k = 1, 
where h is Planck's constant divided by 27r, c is the speed of light, and k is the Boltzmann 
constant. However, for clarity, we reinsert c in the units of quantities whose values are given 
in the text or tables. 

Integration of Eq. (0) over spacetime leads to the Cooper-Frye formula for the direct 
contribution to the invariant one-particle multiplicity distribution namely 

rf^A^dir 1 d^Ndir 2 J + 1 /• 3 p'' 



'^"^^^ dp^ 2iTmt dy drrit {2it)^ is exp{[p ■ v{x) — p{x)]/T{x)} ^ 1 ' ^ ^ 



where E denotes the particle's energy, rrit = \/rn? + p^ its transverse mass, and y its rapidity. 



The quantity m is the particle's rest mass, and pt = \Jpx^ + Py^ is its transverse momentum. 
We assume that the source is boost invariant within the limited region between its two ends 



2^ , |23[| , and that it starts expanding from an infinitesimally thin disk at time t = 0. The 
transverse velocity at any point on the freeze-out hypersurface is assumed to be linear in 
the transverse coordinate p. 



For a particular type of particle, the two-particle correlation function is given by P, |2^j25 



^ '^^ P{pi)P{p2) [/rf%5(a;,pi)][/d%5(x,p2)] ' 

where K = ^{pi + P2) is one-half the pair four-momentum and g = pi — p2 is the pair four- 
momentum difference. The plus sign applies to bosons and the minus sign to fermions, and 
the quantity A specifies the fraction of particles of this type that are produced incoherently. 
The freeze-out hypersurface is specified by 

= l + a,{p/R,f ' 

where Tf is the constant proper time at which freeze-out is assumed to occur along the 
symmetry axis of the source and i?t is the transverse radius of the source at the beginning 
of freeze-out. The transverse freeze-out coefficient specifies the radial behavior of the 
freeze-out and is related to the transverse velocity ft and width Ar in proper time during 
which freeze-out occurs. We obtain for this relationship 



Ar = Tf 



(6) 



under the additional assumption that the exterior matter at 2; = that freezes out first has 
been moving with constant transverse velocity ft from time t = until the beginning of 
freeze-out. 



4 



III. DETAILED ANALYSIS OF THE REACTION Si + Au AT piabM = 14 6 GeV/c 



We now use the nine-parameter expanding source model described in Sec. y to perform a 
detailed analysis of nearly central collisions in the reaction Si + Au at piabM = 14.6 GeV/c, 
for which excellent experimental data were collected in Experiment E-802 ||lO|,|ll| at the 
Alternating Gradient Synchrotron of the Brookhaven National Laboratory. We first describe 
the results of a simultaneous consideration of the pion and kaon one-particle and correlation 
data and then examine various subsets of the data and determine the effect of alternative 
assumptions on our results. 



A. Simultaneous consideration of pion and kaon one-particle and correlation data 

In Experiment E-802, invariant 7r+, vr^, , and K~ one-particle multiplicity distribu- 
tions ll^] and TT"'" and two-particle correlations [|ll[] were measured for the central 7% of 



collisions in the reaction Si + Au at puh/A = 14.6 GeV/c. The nine adjustable parameters 
of our expanding source model have been determined by minimizing with a total of 1416 
data points for the six types of data considered, so the number of degrees of freedom u is 
1407. The error for each point is calculated as the square root of the sum of the squares of its 
statistical error and its systematic error, with a systematic error of 15% for tt"*", tt", and 
one-particle multiplicity distributions, 20% for the K~ one-particle multiplicity distribution, 
and zero for vr"*" and two-particle correlations |10,11|. The resulting value of is 1484.6, 



which corresponds to an acceptable value of x^/v = 1.055. The values of nine independent 
freeze-out properties determined this way, along with their uncertainties at 99% confidence 
limits on all quantities considered jointly, are given in Table |. These values and their un- 
certainties were obtained earlier in Ref . . The quantity hq appearing in Table | and in the 
remaining tables denotes normal nuclear density, whose value is calculated from the nuclear 
radius constant by means of Uq = 3/(47rro'^) = 3/[47r(1.16 fm)^] = 0.153 fm~^. 

TABLE I. Nine independent source freeze-out properties for central collisions of Si + Au 
^-t Piab/^ = 14.6 GeV/c resulting from the simultaneous consideration of the pion and kaon 
one-particle and correlation data. The value used for normal nuclear density no is 0.153 fm^'^. 





Value and uncertainty 


Property 


at 99% confidence 



Central baryon density u/uq 0.145 I0045 

Nuclear temperature T (MeV) 92.9 ± 4.4 

Transverse collective velocity ft (c) 0.683 it 0.048 

Longitudinal collective velocity Vi (c) 0.900 ^0 029 

Source velocity (c) 0.875 to^\l 

Transverse radius Rt (fm) 8.0 it 1.6 

Longitudinal proper time Tf (fm/c) 8.2 it 2.2 

Width in proper time At (fm/c) 5.9 6 

Pion incoherence fraction A^^ 0.65 it 0.11 
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B. Separate consideration of one-particle data and correlation data 



We next use our expanding source model to analyze the invariant 7r+, vr^, , and K~ 
one-particle multiplicity distributions alone, for which there are a total of 656 data points. 
Because the pion incoherence fraction A^r does not enter in the expression for one-particle 
multiplicity distributions, there are only eight parameters in this case, so the number of 
degrees of freedom v is 648. The resulting value of is 692.5, which corresponds to a 
value of x^l^^ = 1-069. The values of eight independent freeze-out properties determined 
this way, along with their uncertainties at 99% confidence limits on all quantities considered 
jointly, are given in the second column of Table It is seen that one-particle multiplicity 
distributions alone determine the freeze-out central baryon density, nuclear temperature, 
transverse collective velocity, longitudinal collective velocity, and source velocity moderately 
well, but contain almost no information concerning the transverse radius, longitudinal proper 
time, and width in proper time. 



TABLE II. Effect on nine independent source freeze-out properties for central collisions of 
Si + Au at piah/A = 14.6 GeV/c of considering separately the one-particle data and the correlation 
data. The symbols denoting the freeze-out properties are defined in Table |. 





Value and uncertainty at 99% confidence 


Property 


One-particle data 


Correlation data 


n/no 


r. ..t- +0.079 
U.i40 _o.057 


'-■^ -1.1 


T (MeV) 


92.9 ± 4.7 


70 


■^t (c) 


0.700 ± 0.095 


0.95 


Vi (c) 


0.904 t'o-'of. 


0.92 10:0« 


^^s (c) 


0.872 t'oTr 


0.96 IHi 


Rt (fm) 


q +20 
^ -9 


9.6 ± 4.9 


Tf (fm/c) 


7 +23 
' -7 


9.6 ± 7.2 


At (fm/c) 


o q +2.4 
"^•y -3.9 


-6.5 






0.75 ± 0.18 



In an analogous study, we use our expanding source model to analyze the vr"*" and two- 
particle correlations alone, for which there are a total of 760 data points, or 751 degrees of 
freedom u. The resulting value of is 756.0, which corresponds to a highly acceptable value 
of x"^/^ = 1.007. The values of nine independent freeze-out properties determined this way, 
along with their uncertainties at 99% confidence limits on all quantities considered jointly, 
are given in the third column of Table 0. It is seen that two-particle correlations alone 
determine the transverse radius, longitudinal proper time, width in proper time, and pion 
incoherence fraction fairly well, but contain almost no information concerning the central 
baryon density, nuclear temperature, transverse collective velocity, longitudinal collective 
velocity, and source velocity. 
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C. Separate consideration of pion data and kaon data 



In our next study, we use our expanding source model to analyze only the data for 
pions, namely the invariant vr"^ and tt~ one-particle multiplicity distributions and the vr"*" 
two-particle correlations. For this case, there are a total of 934 data points, or 925 degrees 
of freedom u. The resulting value of is 959.9, which corresponds to an acceptable value 
of x^/^ = 1.038. The values of nine independent freeze-out properties determined this way, 
along with their uncertainties at 99% confidence limits on all quantities considered jointly. 



are given in the second column of Table |T|. It is seen that when pions alone are considered, 
the freeze-out occurs somewhat later and at a lower temperature than when both pions and 
kaons are considered (Table 1). However, the 99% confidence-level error bars associated with 
the two freeze-outs overlap. 

TABLE III. Effect on nine independent source freeze-out properties for central collisions of 
Si + Au at pia.h/A = 14.6 GeV/c of considering separately the pion data and the kaon data. The 
symbols denoting the freeze-out properties are defined in Table |. 





Value and uncertainty at 9 


9% confidence 


Property 


Pion data 


Kaon data 


n/no 


'^■^^ -0.01 


26 

-0.12 


T (MeV) 


78 ± 20 


101 ± 10 


vt (c) 


0.77 ± 0.14 


0.62 ± 0.14 


ve (c) 


0.920 


0.89 t'ofo 


Vs (c) 


0.880 trolt 


0.847 


Rt (fm) 


9.5 ± 2.7 


7.1 ± 2.5 


Tf (fm/c) 


11.1 ± 5.2 


6.6 ± 3.8 


At (fm/c) 


S A +7-4 
^•4 -4.6 


2 9 +'^-^ 
■^■^ -2.9 




0.75 ± 0.16 





In an analogous study, we use our expanding source model to analyze only the data for 
kaons, namely the invariant and K~ one-particle multiplicity distributions and the 
two-particle correlations, for which there are a total of 482 data points. Because the pion 
incoherence fraction X.^ does not enter when kaons alone are considered, there are only eight 
parameters in this case, so the number of degrees of freedom u is 474. The resulting value of 

is 441.5, which corresponds to a highly acceptable value of x'^/^ = 0.931. The values of 
eight independent freeze-out properties determined this way, along with their uncertainties 
at 99% confidence limits on all quantities considered jointly, are given in the third column of 



Table [ITI[ It is seen that when kaons alone are considered, the freeze-out occurs somewhat 
earlier and at a higher temperature than when both pions and kaons are considered (Table |). 
However, the 99% confidence-level error bars associated with these two freeze-outs, as well 
as with the separate pion freeze-out and kaon freeze-out, overlap. 
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D. Constraining and generalizing the freeze-out hypersurface 



We next consider the effect of constraining the transverse freeze-out to the same source 
time for all points with the same longitudinal position, which is accomplished by holding 
the transverse freeze-out coefficient appearing in Eq. (|^) fixed at zero. We once again 
use all 1416 data points for our six types of pion and kaon one-particle and correlation data. 
Because at is held fixed, there are only eight parameters in this case, so the number of 
degrees of freedom v is 1408. The resulting value of is 1519.4, which corresponds to a 
value of x^/^ = 1.079. The values of nine independent freeze-out properties determined 
this way, along with their uncertainties at 99% confidence limits, are given in the second 



column of Table Because of the strong dependence of the width in proper time Ar 
upon the transverse freeze-out coefficient at, which is held fixed at zero, the value of Ar and 
its uncertainty are anomalously small in this case. We see that constraining the freeze-out 
hypersurface in this way affects the remaining eight quantities only to within their 99% 
confidence-level error bars. 

TABLE IV. Effect on nine independent source freeze-out properties for central collisions of 
Si + Au at piab/^ = 14.6 GeV/c of constraining tfie transverse freeze-out to the same source time 
for all points with the same longitudinal position (transverse freeze-out coefficient at fixed at zero) 
and of allowing a more flexible freeze-out in the longitudinal direction (longitudinal freeze-out 
coefficient fixed at 0.20 c~^). The symbols denoting the freeze-out properties are defined in 
Table |. 





Value and uncertainty at 


99% confidence 


Property 


at fixed at zero 


a£ fixed at 0.20 


n/no 


r. . on +0.058 
U.iOJ _o.042 


r. . +0.061 

-0.043 


T (MeV) 


92.6 ± 4.2 


92.5 ± 4.4 


■^t (c) 


0.721 ± 0.033 


0.663 ± 0.049 


vi (c) 


n nnq +0.022 
u.yuo _o.o27 


925 


Vs (c) 


n 87^^ +0.014 


n ii71 +0.014 
0.8/3 _ooi6 


Rt (fm) 


9.5 ± 1.2 


7.7 ± 1.5 


Tf (fm/c) 


6.3 ± 1.4 


9.7 ± 2.2 


At (fm/c) 


1.94 


7 Q +4.5 
1.0 _2 7 




0.64 ± 0.11 


0.57 ± 0.10 



In a related study, we consider the effect of allowing a more flexible freeze-out in the 
longitudinal direction. To accomplish this, we generalize Eq. (^) to 

where the longitudinal freeze-out coefficient ai allows the freeze-out along the symmetry 
axis of the source to occur with a dependence on longitudinal distance z that is different 
from that corresponding to a constant proper time (so that Tf no longer has this physical 
interpretation). By minimizing with a total of 1416 data points for the six types of data 
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considered, we determined that the minimum in occurs at a value of the longitudinal 
freeze-out coefficient = 0.20 c~^. The resulting freeze-out hypersurface is similar to 
that obtained in nuclear fluid-dynamical calculations by Schlei which was the original 
motivation for this generalization. We then held ae flxed at 0.20 c~^, so that once again 
there are nine parameters and 1407 degrees of freedom. The resulting value of is 1468.5, 
which corresponds to an acceptable value of x^/^ = 1.044. The values of nine independent 
freeze-out properties determined this way, along with their uncertainties at 99% confldence 



limits, are given in the third column of Table IV. It is seen that the primary effect of this 



generalized freeze-out hypersurface is to increase somewhat the values of Tf and At, but 
that the precise shape of the freeze-out hypersurface remains relatively unimportant. 

E. Use of unnormalized one-particle data and inclusion of proton one-particle data 

Because many analyses have been performed with unnormalized one-particle data 



1^-|T8[, we now consider the effect of regarding the pion and kaon one-particle data to 
be unnormalized. Once again, we use all 1416 data points for our six types of pion and kaon 
one-particle and correlation data. Because the pion and kaon one-particle normalization 
constants are allowed to vary, there are 11 parameters in this case, so the number of degrees 
of freedom u is 1405. The resulting value of is 1483.2, which corresponds to a value of 
X^/z/ = 1.056. The values of nine independent freeze-out properties determined this way, 
along with their uncertainties at 99% confldence limits, are given in the second column of 
Table 0. It is seen that the extracted temperature increases slightly, but that its uncertainty 
increases substantially. 

TABLE V. Effect on nine independent source freeze-out properties for central collisions of 
Si + Au at p\ah/A = 14.6 GeV/c of regarding the pion and kaon one-particle data to be unnor- 
malized and of including proton one-particle data (which are contaminated by spectator protons). 
The symbols denoting the freeze-out properties are defined in Table |. 





Value and uncertainty at 99% confidence 




Unnormalized 


Proton one-particle 


Property 


one-particle data 


data included 


n/no 


n 19 

'-'•■'-^ -0.18 


n 140 +0-039 


T (MeV) 


98 ± 27 


97.0 ± 4.0 


vt (c) 


0.65 ± 0.19 


0.614 ± 0.033 


V£ (c) 


n 89S 

u.oyo _o.063 


n Q29 +0-016 
ti.y^y -0.021 


t^s (c) 


n 87f^ +0.016 
U-«'5 _o.018 


U.OOZ _Q 022 


Rt (fm) 


7.8 ± 2.4 


7.4 ± 1.4 


Tf (fm/c) 


8.1 ± 2.7 


8.5 ± 2.0 


At (fm/c) 


6 1 +^-'^ 

"•J^ -3.0 


6.6 tl-^j 




0.66 ± 0.15 


0.66 ± 0.11 



Because proton one-particle data are contaminated by the presence of spectator protons, 
we have not included them in our analysis thus far. However, because they are frequently 
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included in other analyses [^J-|T8[ , we now consider the effect of including 331 data points for 
the proton one-particle multiplicity distribution corresponding to the central 7% of collisions 
in the same reaction that we have been considering |T0[. A systematic error of 15% is used 



also for the proton one-particle multiplicity distribution. With a total of 1747 data points for 
the seven types of data considered and nine adjustable parameters, the number of degrees 
of freedom v is 1738 in this case. The resulting value of is 2546.6, which corresponds to 
an unacceptably large value of jv = 1.465. The probability that a perfect model would 
have resulted in a value of least as large as that found here is the incredibly small 

value 1.1 X 10~^^. Nevertheless, for completeness, we give in the third column of Table [V| 
the values of nine independent freeze-out properties determined this way, along with their 
uncertainties at 99% confidence limits. 

IV. FLAWS IN PREVIOUS ANALYSES 

The results of the above detailed analyses indicate that the freeze-out temperature is less 
than 100 MeV and that both the longitudinal and transverse collective velocities — which are 
anti-correlated with the temperature — are substantial. Similar conclusions concerning a low 
freeze-out temperature have also been reached in Refs. [|T2,r3|. However, other analyses 



l^-p!8[ have yielded a much higher freeze-out temperature of approximately 140 MeV for 



both this reaction and other reactions involving heavier projectiles and/or higher bombard- 
ing energies. In order to reconcile this serious discrepancy, we now examine the features 
in these analyses that erroneously led them to the conclusion of a much higher freeze-out 
temperature. These analyses fall into two major classes, which we consider in turn. 

A. Neglect of relativity in extrapolation of slope parameters to zero particle mass 

One type of analysis ]I^,p!B[] was based upon the extrapolation to zero particle mass of 



extracted slope parameters characterizing the dependence of unnormalized transverse one- 
particle multiplicity distributions upon transverse mass. For a given reaction and type of 
particle, this transverse one-particle multiplicity distribution was represented by the expres- 
sion]^ 

±^.^expf-^l. (8) 



rrit drrit \ 



where A is an arbitrary normalization constant and T^s is the extracted slope parameter. 
Values of Tgfj were extracted in this way for six types of particles originating from three 
separate reactions, namely 7r+, tt", K~^, K~ , p, and p originating from the reaction p + p 
at center-of-mass energy ^/s = 23 GeV, from the 10% most central collisions in the reaction 



2 



To facilitate comparisons with our own expressions, we have transformed the notation used in 



Refs. [14, 15 1 to that used here. 
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S + S at piah/A = 200 GeV/c, and from the 6.4% most central collisions in the reaction 
Pb + Pb at piabM = 158 GeV/c. 

As we will see below, the values of these extracted slope parameters contain valuable 
information, but they were unfortunately analyzed in Refs. [p!4| , [l5| in terms of the erroneous 
equation 



-eff 



T + mf ^ 



(9) 



where T is the nuclear temperature (whose value we are trying to determine) and v is the 
average transverse collective velocity of the expanding matter from which the particle orig- 
inated. Alas, this equation neglects relativity — even though these are relativistic collisions! 
It was introduced on page 182c of Ref. with the phrase "One may empirically guess a 
relationship between the slope parameter and particle mass," whereas the words describing 
the same equation on page 2082 of Ref. |]T3| are "The correlation between the slope parame- 
ter and particle mass m may be described qualitatively by the relationship ..." On the basis 
of the erroneous Eq. (|), the extrapolation in Ref. ||15| of the extracted slope parameters to 
zero particle mass yielded the result T ^ 140 ±15 MeV. 

In the limit in which the particle velocity is large compared to the average collective 
velocity and with the aid of other simplifying assumptions and approximations,^ the cor- 
rect relationship between slope parameter, nuclear temperature, particle mass, and average 
collective velocity can be easily derived from the relativistically correct Eq. (^. With the 
neglect of contributions from resonance decays, the neglect of the =f1 appearing in the 
denominator of Eq. (|^), the assumption of a constant freeze-out temperature, and the as- 
sumption that freeze-out occurs at a constant time t in the source frame, Eq. leads 
to 



1 (fN 
rrit dydrrit 



/ d\ exp 


p ■ f (x) 


= A'E J d^xexpj 


T 


JV 





7(x)[E-p-v(x)] 



(10) 



where A' is a different arbitrary normalization constant from the one appearing in Eq. (^, 
the subscript V on the integral denotes the spatial integration limits for the source, and the 
position-dependent Lorentz factor 7(x) = l/-y/l — v(x)-v(x). 

By introducing an average collective velocity v in the integrations in Eq. ([T0|), taking the 



limit in which the particle velocity is large compared to the collective velocity, specializing 
to the transverse direction, and neglecting the pre-exponential E dependence, we are led to 



1 dN 



rrit drrit 



y4exp 



'J {rrit - Ptv) 
T 



A exp 



rrit — v^m,^ — 



y2 



(11) 



To obtain the relationship between the slope parameter, nuclear temperature, particle mass, 
and average collective velocity, we equate the derivatives with respect to rrit of Eqs. (§) and 
(0), which leads to 



^The expanding source model developed in Refs. |^,^ and used here does not require that the 
particle velocity be large compared to the average collective velocity and does not utilize these 
other simplifying assumptions and approximations. 
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An analogous relationship has also been obtained by Siemens and Rasmussen for the 
case of a blast wave produced by the explosion of a spherically symmetric fireball. 
In the limit of zero particle mass, Eq. (O) reduces to 



T = Tji^. (13) 



which agrees with the result obtained by Schnedermann, Sollfrank, and Heinz |29|j30|] for the 
case of cylindrical symmetry. With a typical value of 0.4 c for the average collective velocity 
V and the limiting value of T^e ~ 140 ± 15 MeV obtained in Ref. [1^ by extrapolating 
slope parameters to zero particle mass, Eq. (|13]) yields T ^ 92 ± 10 MeV for the nuclear 
temperature. 



B. Accumulation of effects from several approximations made in a thermal model 



Another type of analysis [p!^-[18[| utilized the thermal model of Schnedermann, Sollfrank, 



and Heinz [pQjpOl to extract the nuclear temperature and transverse surface collective veloc- 
ity from unnormalized experimental transverse one-particle multiplicity distributions. An 
accumulation of effects from several approximations led to a somewhat higher temperature 
than we have found in our expanding source model PJ^. These approximations include 
the neglect of contributions from resonance decays, the neglect of the =f1 appearing in the 
denominator of Eq. (^), the neglect of the coupling of the transverse motion to the longitu- 
dinal motion, and — most importantly — the neglect of information contained in the absolute 
normalization of the multiplicity distributions. The accumulation of effects from these ap- 
proximations was responsible for the conclusion on page 2083 of Ref. |jl5[ that "Within a 
temperature range 100 < T < 150 MeV, the fits are equally good." Clearly, the use of un- 
normalized experimental transverse one-particle multiplicity distributions in such a thermal 
model cannot be expected to provide a definitive determination of the nuclear temperature 
at freeze-out. 



V. SUMMARY AND CONCLUSION 

We have used a nine-parameter expanding source model that includes special relativity, 
quantum statistics, resonance decays, and freeze-out on a realistic hypersurface in space- 
time to analyze in detail invariant vr"*", vr", K~^, and K~ one-particle multiplicity distribu- 
tions and 7i~^ and two-particle correlations in nearly central collisions of Si + Au at 
PiabM = 14.6 GeV/c. By considering separately the one-particle data and the correlation 
data, we found that the central baryon density, nuclear temperature, transverse collective 
velocity, longitudinal collective velocity, and source velocity are determined primarily by 
one-particle multiplicity distributions and that the transverse radius, longitudinal proper 
time, width in proper time, and pion incoherence fraction are determined primarily by two- 
particle correlations. By considering separately the pion data and the kaon data, we found 
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that although the pion freeze-out occurs somewhat later than the kaon freeze-out, the 99% 
confidence-level error bars associated with the two freeze-outs overlap. 

By constraining the transverse freeze-out to the same source time for all points with 
the same longitudinal position and by allowing a more flexible freeze-out in the longitudi- 
nal direction, we found that the precise shape of the freeze-out hypersurface is relatively 
unimportant. By regarding the pion and kaon one-particle data to be unnormalized, we 
found that the nuclear temperature increases slightly, but that its uncertainty increases 
substantially. By including proton one-particle data (which are contaminated by spectator 
protons), we found that the nuclear temperature increases slightly. These detailed studies 
confirm our earlier conclusion based on the simultaneous consideration of the pion and 
kaon one-particle and correlation data that the freeze-out temperature is less than 100 MeV 
and that both the longitudinal and transverse collective velocities — which are anti-correlated 
with the temperature — are substantial. 

We also discussed the flaws in previous analyses that yielded a much higher freeze-out 
temperature of approximately 140 MeV for both this reaction and other reactions involving 
heavier projectiles and/or higher bombarding energies. One type of analysis was based upon 
the use of an erroneous equation that neglects relativity to extrapolate slope parameters to 
zero particle mass. Another type of analysis utilized a thermal model in which there was an 
accumulation of effects from several approximations. 

The future should witness the arrival of much new data on invariant one-particle multi- 
plicity distributions and two-particle correlations as functions of bombarding energy and/or 
size of the colliding nuclei. The proper analysis of these data in terms of a realistic model 
could yield accurate values for the density, temperature, collective velocity, size, and other 
properties of the expanding matter as it freezes out into a collection of noninteracting 
hadrons. A sharp discontinuity in the value of one or more of these properties could con- 
ceivably be the long-awaited signal for the formation of a quark-gluon plasma or other new 
physics. 
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